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Simulating Dense Matter 
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^ . I review the Sign Problem hindering lattice QCD simulations of dense baryonic matter, 

^ focussing where possible on its physical relevance. The possibility of avoiding the Sign 

' Problem via a duality transformation is also briefly considered. Finally, I review evidence 

for deconfinement at non-zero quark density in recent simulations of Two Color QCD. 
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§1. Motivation 



What is the nature of the QCD ground state in the limit Hb/T 3> 1, where 
T is temperature and \xb the baryon chemical potential? The insight that diquark 
Cooper pair condensation in the color anti-triplet channel is naturally promoted by 
[*"»« ! one-gluon exchange suggests that in this regime QCD is a color superconductor. At 

asymptotic densities \xb — ► oo where weak-coupling methods can be trusted, the 
favoured ground state of QCD with three light quark flavors *' exhibits Color-Flavor 
Locking, with spontaneous symmetry breaking pattern 

O 5Z7(3) C ® SU(3) £ ® SU(3) fl ® U(l) fl ® £7(l)g SU(3) 4 ® U(l) & (IT) 



Italicised and regular fonts denote respectively local and global symmetries of mass- 
less quarks: since both are spontaneously broken the CFL phase is simultaneously 
superconducting and superfluid. 

At the densities available in stellar cores, the QCD coupling g(/xe) is no longer 
small, making reliable calculation difficult. In matter with [ib ~ 0{m s ) pairing 
may only take place between u and d quarks, and further non-trivial constraints are 
imposed by requirements of charge- and color-neutrality. Model approaches have 
predicted many exotic scenarios, such as gap less superconductivity, mixed states of 
normal and superconducting matter, and crystalline LOFF phases. 1 ) The issue of 
which is the true ground state is ideally resolved, of course, by a systematic non- 
perturbative lattice QCD calculation, as suggested by the talk's title. It is worth 
recalling, however, that the most urgent question about quark matter is whether it 
exists at all in our universe inside compact stars, or whether the star would have 
collapsed into a black hole before the required core density can be attained. To settle 
this theoretically we need to solve the Tolman-Oppenheimer-Volkoff equations for 
relativistic stellar structure, which requires quantitative knowledge of the equation 
of state, ie quark density n q , pressure p and energy density e as functions of (is for 
all hb > \iBo-, where fiBo ~ 924MeV is the onset value corresponding to self-bound 
nuclear matter. This issue, surely, is the first goal of lattice QCD with 

*' Astrophysical arguments suggest that QCD matter with four or more quark flavors cannot 
form a stable gravitationally bound system. 
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§2. The Sign Problem, and why we need it 

Let me remind you why this problem was not solved years ago. In Euclidean 
metric the QCD Lagrangian is written 

Cqcd = i^Mif; + -FnuF^u with M(ji) = p[A] + ^70 + m, (2-1) 

where \i = hb/N c is the quark chemical potential. It is straightforward to show 
75M(/x)75 = Mt(— n) implying detM(fj) = (detM(— m))*, so that the path integral 
measure is not positive definite for [i 7^ 0. This is not solely an issue for fermions; 
it can be traced to the explict breaking of time reversal symmetry by the term 
with fj, 7^ 0, which in Euclidean metric corresponds to a breaking of the symmetry 
under i 1— > —i. The consequences are drastic; Monte Carlo importance sampling, the 
mainstay of lattice QCD, becomes ineffective. 

To see why consider the formal solution to the Sign Problem known as reweight- 
ing. Here the phase of detM is treated as an observable, and expectation values 
defined by 

= ((Oarg(detM))) 
[ ' ((arg(detM))) ' 1 ' 

with ((. . .)) defined using the positive measure |detM|e _5i,oson . Unfortunately, both 
numerator and denominator of (2-2) are exponentially suppressed as V — > 00, eg: 

((arg(detM))) = |*™ = exp(-AF) ~ exp(-#V) (2-3) 

6 fake 

where in the last step we assume the free energy F is extensive. On general grounds 
we expect any signal for (O) to be overwhelmed by statistical noise in the thermo- 
dynamic limit. 



Fig. 1. Two polymers of opposite sign. 

It is instructive and surprisingly easy to introduce a Sign Problem into QCD at 
H = 0. Consider the polymer representation 2 ^ for the QCD partition function: 

Z QCD = jDUdetM[U;m]e- Sw[u] 

oc JdU Y,( 2m ) Nm (~ l ) Nr Il(( tr Ylie)(^ u r))e' Sw[U] . (2-4) 



{C} rec eer 

Each non-vanishing term in the expansion of the determinant is represented by a 
partition C of the lattice into N m monomers, dimers and Np polymers, the latter 
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defined as oriented closed paths of links. In the strong-coupling limit where only 
monomers and dimers contribute each term is positive. Polymers contribute not 
only a Wilson loop Up to the effective action but also a shape-dependent sign factor 
— tr]^ gr 7^, where 7^ = ±7^ depending on whether the link £ points along ±p. 
The resulting Sign Problem makes simulation of even the non-interacting system 
difficult. 3 ^ At weak gauge coupling it is tempting to interpret the polymers as quark 
worldlines, but note that the overall sign contributed by C can be changed by the 
innocent-looking flip shown in Fig. 1. It is difficult to believe that in this limit there 
can be any significant correlation between the sign of C and long-range physics. 

The positive measure detM^M used in practical fermion algorithms describes 
color triplet quarks q and color anti-triplet conjugate quarks q c . There are thus 
gauge-invariant qq c bound states with baryon number B 7^ 0. At p = we are con- 
tent to consider these states as extra "mesons" and move on. Once p 7^ 0, however, 
this position is untenable. The lightest baryon in this model's spectrum is degener- 
ate with the pion, so that there is an unphysical onset transition between vacuum 
and baryonic matter at p Q ~ ^m^. Only calculations performed with the correct 
complex measure det 2 M can yield cancellations among configurations with differing 
phases, which nullify the effect of qq c states and postpone the onset transition to the 
phenomenologically-observed p Q ~ ^mjy. For the vacuum to persist as the correct 
ground state in the range p G (^rn^, ^mjv) it seems we actually need a Sign Problem. 

This cancellation has been numerically verified in simulations of Two Color QCD 
(QC2D) with a single staggered quark flavor in the adjoint representation, where it 
was found that the signal for a fake transition at p ~ ^m^, to a superfluid phase 
whose order parameter (qq) vanishes identically due to the Pauli principle, went away 
once the sign of the determinant was correctly taken into account. 4 ) More recently a 
visualisation of the Sign Problem has emerged from an analytic solution of a random 
matrix model with the same global symmetries as QCD, corresponding to the so- 
called mesoscopic limit of V — ► 00 with m^f^V fixed. 5 ) The chiral condensate can 
be expressed in terms of the distribution p(z) of eigenvalues of M — m in the complex 
plane: 

(W>) = lim lim V- 1 [ dzdz p ( Z,m ' ,fi \ (2-5) 

m-»0V-«x> J Z + m 

For p = or the quenched limit Nf = 0, p is real, but in general it is a complex- 
valued function. The explicit solution for p > m n /2 shows that for Re z > m, p 
develops an oscillatory structure, with a characteristic wavelength of 0(V~ l ) and 
amplitude 0(e v ). Any function calculated using a formula such as (2-5) must receive 
wildly fluctuating contributions from different regions of the plane, but remarkably, 
it can be shown that the result behaves entirely in accord with physical expectations, 
namely that {^pip(m)) changes sign at m = 0, but exhibits no sign of discontinuous 
behaviour as p passes though m n /2.^ 

§3. High Density Effective Theory and an Optimistic Conjecture 



For p 2> T, p » Aqcd, QCD is supposed to exist in a deconfined phase of 
degenerate weakly interacting quarks. Is it possible to construct an effective theory 
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\ q = (k + \lk) 

-(p + [ip) 

P 

Fig. 2. Kinematics of quasiquark-gluon scattering at the Fermi surface. 

in terms of quasiparticle degrees of freedom at the Fermi surface, with momenta 
p,k <C /i? Define "fast" (-) and "slow" (+) degrees of freedom via the decomposition: 

ip(x) = exp(i/ix.p) [■;/>+ (x) + ip_(x)] with ip±(p) = -(1 ± a.p)ip(p). (3-1) 

The phase and projection factors ensure that the ip± fields scatter off gluons with 
physical momenta q, according to the kinematics shown in Fig. 2. The quark La- 
grangian in Minkowski metric now reads 7 ) 

V5 +7o (l, f>) v i£> v fl>+ + ^-7o[(l, -Vfibv + M^- + [i>-illD v i> + + h.c], (3-2) 

where A v = e - *' Jxp A i /e*' JX ' p . At tree level the massive ip- field can be integrated 
out to yield 

Chdet = i>+ii\Du^+ + |^ + 7o (7±^) 2 ^+ + C(£> 3 M- 2 ) (3-3) 

with 7| y = ( 7o , P7-P)", 7± = 7" ~ 7j| ■ 

The resulting theory has been used with some success to analyse the color su- 
perconducting phase; however, for our purposes it is more interesting to continue to 
Euclidean metric. 8 ) Since "y\\ u Du is anti-hermitian and satisfies {711^,^^, 75} = 0, it 
is straightforward to show det{^\\ v D v ) is positive definite. There is therefore no Sign 
Problem in the limit /i — > 00! 

I have argued the Sign Problem is intractable almost everywhere in the (fi, T) 
plane as V — > 00 *'. However, it is perhaps possible to distinguish between regions 
such as fi E (^m-K, ^^tin) where we know from Sec. 2 that sign cancellations are both 
subtle and crucial to obtaining physically sensible predictions, from regions where 
the sign fluctuations are not so strongly correlated with long range physics. One 
such region appears to be the upper left-hand corner of the QCD phase diagram 
fi/T ^5 1 where RHIC physics takes place: on finite volumes there is a pleasing 
consistency between approaches based on reweighting (which must inevitably fail in 
the thermodynamic limit) and alternative methods based on analytic continuation. 10 - 1 
While a systematic numerical treatment of HDET is yet to emerge, the previous 
paragraph at least suggests that in the cold dense regime jx/T 3> 1 a "solution" 
of the Sign Problem may not be a crucial component of the physics, and that it 
may be possible to perform controlled calculations on reasonable volumes. Another 

*' An exception appears to be exactly at T — below the fake onset at fi = m^/2. 9 - 1 
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approach, which I will pursue in Sec. 5 below, is to argue that a theory with no Sign 
Problem such as QC2D still models much relevant physics in this regime. 

§4. Are we using the Right Basis? 

Large cancellations between either Feynman diagrams or gauge configurations 
hint at low calculational efficiency. Maybe gauge covariant quarks and gluons are not 
the best degrees of freedom at high density? It is possible in some cases to effect a 
transformation to another basis in which the Sign Problem is absent. An intriguing 
example illustrating this comes from 3d scalar QED, 11 ) with Lagrangian 

Csqed = \f 2 + |Z^| 2 + m 2 \4>\ 2 + A|0| 4 - \e tjk H t F jk , (4-1) 

where H is a real vector source coupled to a generalised B-field. The model has a 
transition separating Coulomb and Higgs phases which is second order for sufficiently 
large A/e 2 . There is a conjectured duality at this critical point with a complex scalar 
field theory described by 

Csft = [{d - eH) k 4?][(d + eH) k 4>] + m 2 \j>\ 2 + A|0| 4 + • • • (4-2) 

The point is that eH^ with e = 2ir/e is a real chemical potential associated with the 
conserved charge density 2Im(</>* d^cp) . The Lagrangian (4-2) is in general complex, 
describing planar Bose-Einstein condensation (BEC) of ^-quanta for eH^ s=s rh. In 
principle, however, it could be studied via simulations of the real action (4-1). 

A more recent example due to Endres 12 ** exploits an exact duality between com- 
plex scalar field theory in d dimensions and a loop gas: 

Z(p) = /iWexp[- £( 2 ^ " <t>>^'°<t> x+ v - ^e-^VJ 

-m 2 Y,€<t>x-Y, V ^M_ 

X X 

oc £ r p Dpe^,o Jj[ e -(2*f2+m"K-v«) Y[l ixA 2p xPx+u )] , (4-3) 

{£} J ° x v 

where p x = \<f) x \, and we have introduced integer- valued link variables l x ^ v governed 
by the constraint d~ ' i v = 0. Once again Z is recast in terms of a functional integral 
over real variables, and the Sign Problem averted. Remarkably, it has proved pos- 
sible to simulate the action (4-3) efficiently, yielding non-trivial results of physical 
interest. 12 ^ 

§5. Two Color Matters 

QCD with gauge group SU(2) and an even number of fundamental quark flavors 
has a real functional measure even once p, 7^ 0, and remains the only dense matter 
system with long-range fundamental interactions amenable to study with orthodox 
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lattice methods. Since q and q live in equivalent representations of the color group, 
hadron multiplets contain both qq mesons and qq, qq baryons. Near the chiral limit 
at // = we expect spontaneous chiral symmetry breaking implying m n <C m p , 
where p denotes the next lightest hadron. The theory's /^-dependence for [i < m p 
can be analysed using chiral perturbation theory 13 ) (xPT). The key result is that 
a second order onset transition occurs at fi Q = \m 1T to a phase with quark charge 
density n q > 0. For Nf = 2 the matter which forms is composed of tightly-bound 
scalar diquarks, which Bose condense to form a gauge-invariant superfluid BEC 
(qq) = (ip tr C^a2T2ip) ^ 0, where the matrices act on spinor, flavor and color indices 
respectively. Since n q — > as fj, \ fi the matter can be arbitrarily dilute. 

The xPT predictions for ($ip(/j,)) and (qq(fJ>)) have been confirmed by lattice 
simulations with staggered fermions. 4 )' 14 ) For our purposes the most relevant pre- 
diction is the equation of state for T = 0, fi > fi Q : 

n xPT = 8N f f^(l-^y (5-1) 
leading to the pressure p x pt = J^ o n q d\i and energy density £ x pt = ~P + pn q : 

p xPT = AN f fl (v 2 + ^- ; txPT = Wffi L 2 - 3^ + 2^) • (5-2) 

This is to be contrasted with another paradigm for cold dense matter, namely a 
degenerate system of weakly interacting quarks populating a Fermi sphere up to 
some maximum momentum kp ~ Ep = fj,, obeying Stefan-Boltzmann (SB) scaling: 

N f N c „ N f N c 4 , , 

nsB = 3?r2 H ; ess = 3psp, = ^ 2 V ■ (5-3) 

The appearance of N c underlines that (5-3) describes a deconfined phase. Superflu- 
idity in this scenario arises from condensation of diquark Cooper pairs within a layer 
of thickness A centred on the Fermi surface, implying (qq) oc Afi 2 . Fig. 3 shows the 




Fig. 3. Comparison of x?T with free quarks. Fig. 4. Equation of state from lattice QC2D. 

ratio of xPT to SB predictions as functions of [ij [i for the choice f 2 = N c /Qir 2 . By 
equating pressures, we naively predict a first-order deconfining transition from BEC 
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to quark matter at Hd/y> ~ 2.3 *'. 

To test whether this prediction is robust we have performed lattice simulations 
of QC2D using Nf = 2 Wilson quarks. 15 ^ Initial runs have been on a 8 3 x 16 lattice 
with lattice spacing (defined via the string tension) a ~ 0.22fm, corresponding to 
T ~ 60MeV. We have not attempted to get particularly close to the chiral limit; 
m^a = 0.79(1) and m^jm, p = 0.80(1). The code's only novelty is the inclusion of a 
diquark source term of the form j(qq + qq), with ja = 0.04 for the most part, which 
both ensures ergodicity and regularises IR fluctuations in the superfluid phase. Fig. 4 
shows the resulting curves for n q , p and e q in the same format as Fig. 3 (open symbols 
denote the j — > extrapolation). There appears to be a transition from confined 
bosonic "nuclear matter" to deconfined quark matter at /i^a ~ 0.65. For large fi, 
n q /nsB ~ p/psb ~ 2, consistent with a bound system having Ep < kp. The claim 
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Fig. 5. Evidence for deconfiricment. 

is supported by Fig. 5, which shows the Polyakov line rise from zero at /id, coincident 
with the superfluid order parameter assuming the scaling (qq) oc fj? consistent with 
Cooper pairing at a Fermi surface. A condensed matter physicist would refer to 
this as a BEC/BCS crossover. Similar conclusions have been reached in a study of 
topological charge susceptibility using staggered fermions with Nf = 8. 16 ^ 
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Fig. 6. Electric gluon De in the static limit. Fig. 7. Magnetic gluon Dm in the static limit. 

A major motivation for studying QC2D is to understand how deconfined quarks 
affect gluodynamics; as argued above, it is at least plausible that the lessons learned 

*^ The apparent transition at [id/ Ho ~ 1-4 can be eliminated by introducing a bag constant to 
stabilise the confined phase at low density. 
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may apply to physical QCD. In any medium with a preferred rest frame, the gluon 
propagator can be decomposed as follows: 

D^(q) = Vff u D M (ql q.q) + V^D E (q 2 , q.q) + V^D L (q 2 , q.q) (5-4) 

with 

t)M _ e _ 1'ilj -pM _ <pM _ n . -p_B _ r _ 1ji.1v _ <r>M . *r,L _ QjjUv 

J J q.q ^ lo+1-l Qo+1-1 

We have used Landau gauge in which = 0. Figs. 6 and 7 plot -De and Z?m in the 
static limit go = as functions of |<f|, for various [i on either side of the deconfinement 
transition. In the electric sector for fia > 0.9 Fig. 6 shows evidence for some Debye 
screening as \q\ — > 0. Deconfinement has a much more dramatic effect in the magnetic 
sector shown in Fig. 7, where in the same limit the propagator is screened by O(50%). 
This is significant because in perturbation theory magnetic gluons are not screened 
in the static limit; indeed, this is at the origin of the celebrated scaling A oc e~#/ 9 
of the color superconducting gap predicted by weak-coupling methods. 17 ) 

In summary, even models with no Sign Problem may hold surprises for us at 
large fi. It would, of course, be nice to compare lattice results for QC2D with other 
non-perturbative approaches such as Schwinger-Dyson equations. My overall feeling, 
though, is that a radical reformulation of non-perturbative QCD is needed before 
numerical approaches can make further headway. 
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